The Physics of CMBR Anisotropies 
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The observed structures in the universe are thought to have arisen from gravitational instability acting on small 
fluctuations generated in the early universe. These spatial fluctuations are imprinted on the CMBR as angular 
anisotropies. The physics which connects initial fluctuations in the early universe to the observed anisotropies is 
fairly well understood, since for most part it involves linear perturbation theory. This makes CMBR anisotropies 
one of the cleanest probes of the initial fluctuations, various cosmological parameters governing their evolution 
and also the geometry of the universe. We review here in a fairly pedagogical manner the physics of the CMBR 
anisotropies and explain the role they play in probing cosmological parameters, especially in the light of the 
latest observations from the WMAP satellite. 



I. INTRODUCTION 



The cosmic microwave background radiation (CMBR) is of fundamental importance in cosmology. Its serendipitous discovery 
by Penzias and Wilson 1 1], gave the first clear indication of an early hot 'Big bang" stage on the evolution of the universe. The 
subsequent verification by host of experiments, culminating in the results of the COBE satellite confirmed that its spectrum is 
very accurately Planckian |2], with a temperature T = 2.725. This is the firmest evidence that the universe was in thermal 
equilibrium at some early stage. Indeed the observed limits on the spectral distortions severely constrain any significant energy 
input into the CMBR below z < 10 7 or so 0|. 

Shortly after its discovery, it was also predicted that the CMBR should show angular variations in its temperature, due to 
photons propagating in an inhomogeneous universe 0|. In the standard picture, the baryonic matter in the early universe was in 
a highly ionized form with radiation strongly coupled to the baryons. As the universe expanded, the matter cooled and atoms 
f^S i formed below about 3000° K. After this epoch the photon mean free path increased to greater than the present Hubble radius, and 
they could free stream to us. These are the photons that we detect in the CMB. They carry information both about the conditions 
at the epoch of their last scattering, as well as processes which affect their propagation subsequently. Fluctuations in the early 
universe result in inhomogeneities on the 'last scattering surface' (henceforth LSS). These inhomogeneities should be seen today 
as angular anisotropies in the temperature of the CMB. Further, the CMB photons are influenced by a number of gravitational 
and scattering effects during their passage from the LSS to the observer. These are also expected to generate additional CMBR 
anisotropies. 

These CMBR anisotropies took a long time to be discovered and its absence in the early observations were beginning to prove 
embarrassing, for theories of structure formation. It was not until 1992 that the temperature anisotropies in the CMBR were 
detected, on large angular scales, by the Differential Microwave Radiometers (DMR) experiment on the COBE satellite |5]. The 
fractional temperature anisotropies are at the level of 10~ 5 and ruled out some of the earlier baryon dominated models, and hot 
J> , dark matter dominated models, but were quite consistent with expectations from latter Cold Dark matter models of structure 
formation |6,0]. 

Since the COBE discovery a large number of expreiments have subsequently probed the CMBR angular anisotropies over a 
variety of angular scales, from degrees to arc minutes (cf. |8] for a recent review). This has culminated in the release of the 
first year all-sky data from the Wilkinson Microwave Anisotropy Probe (WMAP) satellite |9]. These observations, especially 
the 'acoustic oscillations' which are inferred from the anisotropy power spectrum, have led to the confirmation of a popular 
'standard' picture for structure formation; one where an early epoch of inflation generated adiabatic perturbations in a spatially 
flat universe. The observed anisotropy patterns also allow cosmological parameters to be probed with considerable precision, 
especially when combined with other data sets related to the observed inhomogeneous universe II (111 211. It has therefore become 
imperative for the modern cosmologist to understand the physics behind CMBR anisotropies. We review here in a pedagogical 
fashion, the relevant physics of the temperature anisotropies and also briefly mention the polarization of the CMBR. There are 
a large number of reviews lo[l4fl5lflal . and text books 1171 [l8ll on this subject. The authors aim is to present some of these 
ideas in a manner in which he, as a non expert, understood the subject, which may be of use to some! 



II. THE CMB OBSERVABLES 



The CMB is described by its brightness (or intensity) distribution. Since the spectrum of the CMB brightness, seen along any 
direction on the sky n, is very close to thermal, it suffices in most cases to give the temperature T(n). The temperature is very 
nearly uniform with fluctuations AT(n) at the level of 10~ 5 T, after removing a dipole contribution. It is convenient to expand 
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the temperature anisotropies AT(n)/T = 9(n) at the observer in spherical harmonics 
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<3(n) = — (6,0) =Y,ai m Y lm (6,0) (1) 
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with a'[ m — (— \) m ai^ m , since the temperature is a real quantity. 

In the standard picture, the universe is assumed to have evolved from density fluctuations initially described by a Gaussian 
random field, and one can then take 9 to be a Gaussian random field. In this case a; m 's are also Gaussian random variables with 
zero mean and a variance completely described by their power spectrum, 

(a>lm a l'm') — Ci5u'S mm '. (2) 

Here we have assumed also the statistical isotropy of O(n) field because of which the power spectrum is independent of m. 
Theoretical predictions of CMBR anisotropy are then compared with observations by computing the C;'s or the correlation 
function C(a) = (O(n)O(m)), where if we have statistical isotropy, C depends only on cos a = n ■ m. From Eq. (|2ji and the 
addition theorem for the spherical harmonics, we have 

o; i i 

C ( a ) = Y.Y.i^'^YimY^, =J2Ci-^-Pi(coaa). (3) 

Im I'm' I 

The mean-square temperature anisotropy, ((AT) 2 ) = T 2 C(0) is 

((AT) 2 ) = T 2 ^ C; ~ T 2 f + 2 l )Q dlnl (4) 

with the last approximate equality valid for large I, and so I (I + 1 )Ci /2ir is a measure of the power in the temperature anisotropies, 
per logarithmic interval in / space. (We will see below that scale invariant potential perturbations generate anisotropies, due to 
the Sachs-Wolfe effect @|, with a constant 1(1 + 1)C/, which provides one motivation for this particular combination). 

Note that the CMB brightness and hence 9 is also a function of the space-time location (x, 77) of the observer. Here x is the 
conformal spatial and 77 the conformal time co-ordinates respectively (see below). One computes the correlation function C(a) 
predicted by a given theory by taking the ensemble average of (9(xo, ?7o, rl )9(xo, 770, m)). For the statistically isotropic case 
this again only depends on cos a = n • m. Further the Fourier component of 9, for every k mode, often depends on n only 
through k ■ n = \i, where k = k/\k\. One can then conveniently expand 9 in a Fourier, Legendre series as 

/r/ 3 k . 
__ e4k .xo £(_i)«(2i + l) ai (fc,77o)P ; (fc ' n) (5) 

For a homogeneous, isotropic, Gaussian random 9 field, (a;(k, 770)0^, (p, 770)) = (\ai(k, ?7o)| 2 )<5z.z' (27r) 3 <5 3 (k — p), where the 
power spectrum (\ai(k, 770 ) | 2 ) depends only on k = \k\. One then gets 



^W^E 1 / ^fc 3 (h(Mo)| 2 )^iMcosa), (6) 
^ 7T J k Att 



where we have used the addition theorem 

A-jr 

Pl(k-n) = ^ Yi m (n)Y lm (fe) (7) 

m 

and carried out the angular part of the integral over d 3 k using the orthogonality of the Spherical harmonics. Comparing Eq. 
with Eq. (0 we see that 

Q = IJ ^fc 3 (h(fc,%)| 2 ) (8) 

We will use this equation below to calculate Cfs for various cases. One can roughly set up a correspondence between angular 
scale at the observer a, the corresponding I value it refers to in the multipole expansion of (9 2 ) and also the corresponding 
co-moving wavenumber k. One has (a/l°) « (100//) and / w kR* where R* is the comoving angular diameter distance to the 
LSS and is ~ lOh^ 1 Gpc, for a standard ACDM cosmology (see below). 

We show in Figure 1 a plot of the temperature anisotropy AT = TWUj + l)Ci/2n and polarization anisotropy versus /, for 
a standard ACDM cosmology, got by running the publicly available code CMB FAST 1 19]. One sees a number of features in 
such a plot, a flat plateau at low / rising to several peaks and dips, as well as a cut-off at high /. Our aim will be to develop a 
physical understanding of the various features that the figure displays. We now turn to the formalism for computing the C; 's for 
any theory. 
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III. THE BOLTZMANN EQUATION 



The distribution function for photons, f{x l 1 pp) = /(x, p, 77), is defined by giving their number in an infinitesimal phase 
space volume, dN — f(x a ,pp, rj)d 3 x a d 3 pp. Note that we will use Greek letters for purely spatial co-ordinates and Latin ones 
for space-time co-ordinates. We can write dN in a co-ordinate independent way, 

dN= I ^Lfix^p^ptdX^Slp^™} (9) 

J pa V 9 

where g is the determinant of the metric, dY^ — y/—g eijki [dx l A dx^ A dx k ]/3\ is an infinitesimal spacelike hypersurface, pi 
the photon 4-momentum and the delta function ensures that pi is a null vector. To get the simple expression for dN, one chooses 
a time slicing with dE; = (y/— gdx 1 dx 2 dx 3 , 0, 0, 0) and carries out the integral over po retaining only the positive energy part 
of the delta function. This shows explicitly that / is a co-ordinate invariant scalar field. Further, in the absence of collisions, 
both dN and the phase-space volume d 3 x a d 3 pp would be conserved along the photon trajectory and hence also the phase space 
density /. If A is an affine parameter along the null geodesic, we will then have df/dX = 0. On the other hand when collisions 
are present the distribution function will change. This situation is generally handled by introducing a collision term on the RHS 
of the CBE, that is writing df /dX — c(f). Further it is generally convenient to use the time co-ordinate itself, say 77, as the 
independent parameter along the photon trajectory and write df/dX — (dr)/dX)(df Jdrj) = c(/) = (dr]/dX)c(f), One then has 

df_ _ df_ + dx^_df_ + dpad£_ _ 
di] dr/ drj dx a drj dp a 

We look at this equation in the spatially flat, perturbed FRW universe. Its metric in the conformal-Newton gauge is 

ds 2 = a 2 {r]) [(1 + 2(j))dr} 2 - (1 - 2<j)){dx 2 + dy 2 + dz 2 )] (11) 

Here 0(77) is the expansion factor and 77 is the conformal time, related to the proper time by a{r])drj = dt. (We adopt c = 1 
units). We have assumed that a single potential <fi describes the scalar perturbation, which holds if the source of the perturbations 
is a perfect fluid with no off-diagonal components to the energy-momentum tensor. 

Since the photon 4-momentum p l is a null vector we have g lk PiPk — 0. We choose the photon 4-momentum to have compo- 
nents pi = (y/ g aa I g 00 p, —pn a ) where, we have defined the magnitude of the spatial component of (co-variant) momentum, 
p = Y. a p a p a . Also n is a unit vector in the direction of the photon S-momentump". 

Then to linear order in the perturbed potential, 

The geodesic equation for the photons to the linear order in the perturbations is 

^-2 P ^. (13) 



dr/ dx l 



The Boltzmann equation then becomes 



Tr, = ^ + (1 + 20)n ^ + 2p d^dp~ a = c(/) (14) 

An observer at rest in the perturbed FRW universe, has a 4-velocity u l = {l/^/goo, 0). So the energy of the photon measured 
by such an observer is E = piU 1 = p(\ + <p)/a. In the unperturbed FRW universe, the energy simply redshifts with expansion 
with E = p/a. The distribution function for the photons, in the absence of perturbations is then described by the Planck law, 



A (f) 



1 (15) 



exp(p/T) — 1 



Defining the perturbed phase space density fi(x,p, n, 77) = f(x,p,n,rf) ~ fb(p), to linear order the Boltzmann equation 
becomes 

^l+n-Vf 1 -2pn-Vcb^ = c(f) (16) 
where we have replaced / by /& in the term last on the LHS of Eq. H41 and used (dfb/dp a ) = (p a /p)(dfb/dp) = 
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In the perturbed FRW universe we note that both the perturbed trajectory and the effect of collisions (under the Thomson 
scattering approximation) do not depend on the photon energy. This motivates us to define the perturbed phase space density in 
terms of a purely temperature perturbation ST(x, T), n) in /;,. We take 



P 



T + 8T{x,ri,n) 



(17) 



(Note that in such an approximation, we are also neglecting the effects of any spectral distortion). To the first order in ST/T 
one can expand / in Eq. (II 7i to get /j = —p(dfb/dp)(6T/T). Again because both the perturbed trajectory and the effect of 
collisions do not depend on p, we will usually integrate overp's. It is then useful to deal with not the full phase space density, 
but just its associated fractional brightness perturbation defined by 

To appreciate better the meaning of (5T/T) = i/4 let us look at the energy momentum tensor of the photons. This is 

d 4 Pk 



(18) 



rp% 

3 ~ 



p l p 3 2 5{p m p m )f = 



d 3 p a tfPj 
,0 J 



-g J v-9 p 

The energy density in the perturbed FRW universe, with metric given by Eq. dill is 

(1+40) 



p = T o = (l±_^) J p s dpdnf = pR(1 + 4(f)) \l + Ji™ 



(19) 



(20) 



where pr — (An/a 4 ) J p 3 dpfb(p) is the radiation energy density in the absence of perturbations. Let us define io = f i(dQ,/Aii) 
and (5T/T)o — J (ST/T) (df2/47r), the zeroth moments of the perturbed brightness i and temperature (ST/T) respectively, over 
the directions of the photon momenta. The fractional perturbation to the radiation energy density is given by 



P~ PR 
PR 



«0 



ST 
~T 



(21) 



Note that in the conformal Newton gauge the radiation density perturbation has an additional contribution from the perturbed 
potential itself (over and above the contribution from the perturbed distribution function). One may feel that this differs from 
the naive expectation that the radiation energy density to vary as p oc T 4 and hence the 'physical' temperature perturbation go 
as just (l/A)8p/p. Since the energy of a photon seen by an observer at rest in the perturbed FRW universe is E = p(l + <fr)/a 
(see above), one can write the phase space density in the perturbed FRW universe, to linear order as 



/ = fb 



P 



T + ST 



= fb 



Ea 



= fb 



Ea 



(22) 



(T + 8T){1 + <t>) J \T+(5T + q 

This shows that = ST + is indeed the 'physical' temperature perturbation measured by an observer at rest in the perturbed 
FRW universe and that Qq = j ®(d£l/4ir) = Sr/A, as expected. 



IV. THE COLLISION TERM AND THE EQUATION FOR THE PERTURBED BRIGHTNESS 



We now consider the effect of collisions. The process that we wish to take into account is the scattering between photons 
and electrons. In fact to linear order it is sufficient to consider the Thomson scattering limit of negligible energy transfer in the 
electron rest frame. Since the distribution function is a scalar, the effects of collision are most simply calculated by going to 
the electron (or fluid) rest frame and transferring the results to any other co-ordinate frame. Suppose the distribution function 
changes by df = df = drc in the fluid rest frame. Henceforth quantities with an 'overbar' will represent variables in the 
fluid rest frame. Then one can write df /drj — {d,T/drj)c. The differential cross section for Thomson scattering of unpolarized 
radiation is given by 



da 

dri* 



<7y 
47T 



1 



Pa(n-n') 



(23) 



where ctt is the Thomson cross section and n and n' are the unit vectors specifying the direction of the initial and the scattered 
photon momenta in the fluid rest frame. The collision term df/dr = df/dr = c will have a source due to the photons scattered 
into the beam from a direction n' and a sink due to scattering out of the beam. So we have 



c(f) = n e a T j — 



1 



P 2 (n-n') 



[/(?,*»') -/(P,n)] 



(24) 
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where the integration over dil' is over the directions of n'. 

In order to derive the equation satisfied by the brightness perturbation we multiply the linearized Boltzmann equation d!6> by 
p 3 and integrate over p to get 



PR,a 
An 



— + n.Vz + 8n.V0 

or] 



(25) 



We simplify the collision term on the RHS of Eq. i25\ . in Appendix lAl From Eq. i25\ . Eq. (IA4> and Eq. iA5\ the equation 
satisfied by i, to the leading order of the perturbations, is given by 



tt- +n.V« 

orj 



in.Vcj) = n e aT(i 



+ 4n.v+^^F 2m ( 



nk 2r , 



(26) 



The effects of Thomson scattering is to drive the photon distribution such that the RHS of Eq. (I26> would vanish. If the scattering 
cross section had been isotropic, then i would have been driven to io in the fluid rest frame; but in the frame where the fluid 
moves there is a doppler shift. In addition the anisotropy of Thomson scattering introduces the dependence on the quadrupole 
moment of the brightness. The perturbed brightness equation will be used to derive the equation for the CMBR anisotropics, 
and also the dynamics for the baryon photon fluid. 



V. INTEGRAL SOLUTION FOR CMBR ANISOTROPIES 



Consider the perturbed brightness equation Eq. (I26> in Fourier space, in terms of the Fourier coefficients 0(fc, r], n) of the 
'physical' temperature perturbation = ST/T + <j> = i/A + <j), that is 



Q(x,T],n) 



d 3 k 



°e(k,r),n) 



(27) 



Henceforth we shall denote the Fourier transform of any quantity A by A, except for the velocity (v) and potential (<fr) whose 
Fourier transforms are denoted by V and <£> respectively. We also assume that these Fourier co-efficients depend on n only 
through ]i = k ■ n, as will obtain for example for scalar perturbations. In this case one can choose an axis for each k mode such 
that i 2m = for m ^ and F 2m (n)i 2 m/10 = -0 2 -P2(m)/ 2 , where 



e(fe,77, M )=^;(-i)'(2Z + l)e / (fe,7 ? )P / (/i) 



(28) 



is the expansion of in a Legendre series. For scalar perturbations, V is also parallel to the k vector, and son-V = (n-k)V = 
V ]i. Further, it is much more convenient to work with the equation for the combination + $. From the Fourier transform of 
Eq. ( I26> we then have 



(0 + $) + (ik]i + n e cr T a)[Q + $] = n e a T a S{k, jj, fx) + 2$, 



(29) 



where, henceforth an over dot will denote a partial derivative, i.e. / = df /drj, for any /. We have also defined the source 
function 



S(k,r),n) 



©o + $ + Vn- 



p 2 (m)0 2 



(30) 



Suppose we define the differential optical depth to electron scattering dr = n e aTa dr] — n e aTdt. Then one can solve Eq. j29t 
formally as 



©(fc^c/iR^Mo) = [©(%) + <S>(r ll )}e-^") e -r(v , m) 



drj e 



rS(k, r/,n)+2$ 



where we have defined the optical depth to electron scattering between epochs 77 to the present 770 



T iVo,v)= dr]'n e (r]')a T a(r]'), 

Jr) 



(31) 



(32) 
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and t(rf) — dr(r), rj')/dr) = n e (r])(TTQ-{v)- The first term on the RHS of Eq. (13 11 can be neglected by taking a small enough 
initial time r^, because the exponential damping for large optical depths. One can then simplify Eq. ( 13 U to get at the present 
epoch ijo 



We have defined above the visibility function 



d V S(k, ri, n)g{r) Q , v ) e - lk ^°-v) + 2 / drj $ e 

Vi Jvi 



-t(vo,v) — 



n e (r})a T a(r}) exp 



dr]'n e (r)')a T a(ri'} 



(33) 



(34) 



such that .9(770, r))di] gives for every 7/0 the probability that the last scattering of a photon occurred in the interval (r/, r\ + drj). 
Suppose r/o is the conformal time at the present epoch. Then as 77 decreases from 770, the optical depth to electron scattering will 
increase and so will g. However far back into the past when r> l,j will be exponentially damped. So the visibility function 
generally increases as one goes into the past attains a maximum at an 'epoch of last scattering' and decreases exponentially 
thereafter. Its exact behavior of course depends on the evolution of the free electron number density during the recombination 
epoch and also on the subsequent ionization history of the universe. If the universe went through a standard recombination epoch 
with no significant reionization thereafter, then the 'surface of last scattering' is centered at z « 1100 with a very small half 
width Az ss 100. If on the other hand the universe got significantly reionized at high redshifts, as it seems to be indicated by the 
WMAP observations, some fraction r r i ~ 0.17, of the photons will suffer last scattering surface at later epochs. 

We can calculate ai (k, 770) by taking the Legendre transform of both sides of Eq. ( 1331 . Note that the term $(k, 770) on the LHS 
of Eq. (I33> . does not depend on the photon direction and so does not contribute to CMBR anisotropy at all. Using the expansion 
of plane-waves in terms of spherical waves, 



e- ih » x =£(-i)'(2J + l)il(s)fl(M). 



(35) 



and writing /zexp(— i/j,kx) = id(exp(— ifikx))/d(kx), we get 



no 

ai(k,r] Q ) = I drj 5(770,77) 
'0 



drj e 



-T(lO,l)oi 



2$J;(fcA77) 



(36) 



Here ji(kx) is the spherical Bessel function, and j[ denotes a derivative with respect to the argument, A77 = 770 — 77. 

Let us interpret the various terms in Eq. (I36> . This equation shows that anisotropies in the CMBR result from a combination 
of radiation energy density perturbations 0o and potential perturbations $ (the monopole term), on the last scattering surface 
and the doppler effect due to the line of sight component of the baryon velocity V (the dipole term). The anisotropy of the 
Thomson scattering cross section also leads to a dependence on the radiation quadrupole 62. The spherical Bessel function and 
its derivatives in front of these terms project variations in space, at the conformal time 77 around last scattering, to the angular (or 
I) anisotropies at the present epoch 770. (A popular jargon is to say that the monopole, dipole and quadrupole at last scattering 
free stream to produce the higher order multipoles today). These spherical Bessel functions generally peak around kAr/ « I. 
The multipoles I are then probing generally spatial scales with wavenumber k ~ Z/A77 at around last scattering. The visibility 
function weighs the contribution at any conformal time 77 by the probability of last scattering from that epoch. Finally, the last 
term ($ term) shows that any variation of the potential along the line of sight will also lead to CMBR anisotropies, and is usually 
referred to as the integrated Sachs-Wolfe (ISW) effect. 

In the limit of a very narrow LSS at 77 = 77*, For angular scales (and I's) such that, their associated spatial scales at the LSS 
are much larger than the LSS thickness, one can take the variation of the j',s with 77 to be much slower than that of the visibility 
function. In such a narrow LSS approximation, we get 



a ; (fc,77o) = (Go + $W)MkR*) + iV^jKkR*) 



e 2 (rf) (3j'/(kR*) + jtikR*)) 



d^nikAr,), (37) 



where R* = 770 — 77* is the co-moving angular diameter distance to the LSS. Note that due to the presence of e~ T in the in the 
last term, the range of integration is restricted to be from about 77* to the present. The presence of a finite width of the LSS causes 
a contribution to the ISW effect from epochs just around last scattering as well, usually referred to as the early ISW effect. Once 
we calculate the photon brightness, and the baryon velocities at the epochs corresponding to last scattering, one can calculate a/ 
using Eq. ( I36l l and C; from Eq. l|8}. Before considering the dynamics of the the baryon-photon fluid in detail, let us first use the 
Eq. d37l > to calculate the CMBR anisotropies at large angular scales. 
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VI. SACHS-WOLFE EFFECT AND LARGE ANGLE ANISOTROPIES 



We wish to calculate the anisotropies generated at large angular scales (or small values of I), large enough such that the 
associated spatial scales are larger than the Hubble radius at the LSS (i.e. kif <C 1). For such scales one can neglect the 
thickness of the LSS and calculate a; using Eq. d37i . Let us also assume that the universe is spatially flat and that it is matter 
dominated by the time r\ = 77* . The evolution of the gravitational potential <\> is considered in detail in the review by Mukhanov 
et al. 1 20] for a variety of initial conditions, and in various epochs. We will draw upon several of their results below. 



A. Adiabatic perturbations 

Consider first adiabatic (or isentropic) perturbations, for which 5p n /(p n +p n ) is the same initially for all components. (Here 
p n is the pressure of component ri). This condition is preserved by the evolution on super Hubble scales I20B . As we show 
below, it is also preserved in the evolution of the tightly-coupled baryon-photon fluid. For a flat matter dominated universe, the 
potential evolves as + (6/77)0 = 0, which implies that is constant in time, ignoring the decaying mode 1 20] . A detailed 
calculation starting from an initial potential perturbation 0$ and following the evolution of the adiabatic mode from radiation 
era through the matter radiation equality gives, = (9/10)0,; = 0o- The perturbed Einstein equation also gives for the matter 
density perturbation 1 20], 5 m = —20 + (rj 1 /6)V 2 0, and v = — (1/3)?7V0. For adiabatic perturbations 5r = (4/3)<5 m . So 
©o = <W 4 = <W3 = -(2/3)0 + (t7 2 /18)V 2 — > -(2/3)0, for large scales, such that kr\ < 1. So the Fourier co-efficient 
0o + $ = —(2/3)$ + $ — <I>/3. Further since V = — (i/S)(krf)$>, this dipole term in Eq. d37l > is negligible compared to 
the monopole term 0o + 3>. Also because of tight coupling and negligible thickness to the LSS there is negligible quadrupole 
component to 02 for k-q <C 1. On super Hubble scales, for adiabatic perturbations one then has 

., UD ~ nSW _2 [ dkk 3 (\<S> (k)\< 



ai = a? w = -#o MkR*)\ C l =-] T 9 (38) 

The above C/, which describes the CMBR anisotropies on large scales due to initially adiabatic potential perturbations, was first 
derived by Sachs and Wolfe |4] and is referred to as the Sachs-Wolfe effect. For a power law spectrum of potential perturbations, 
with A 2 = fc 3 |$ (fc)| 2 /(2^ 2 ) = A\ (fcoXfc/fco)"'- 1 , one gets 

lSW _ 2A l ( 1 A 2"^ 4 r(3 - n s )T((2l + n s - l)/2) 



1 9 \k R*J r 2 ((4-n 5 )/2)r((2/ + 5-n s )/2)- k ; 

(In the above equation T(x) is the usual gamma function). In theories of inflation, one obtains a nearly scale invariant spectrum 
corresponding to n s = 1. For this case, one gets a constant 



Z(Z + l)Cf w (A 



2?r V 3 



(40) 



It is this constancy of (1(1 + l)Ci)/2w for scale invariant spectra that motivates workers in the field to use this combination to 
present their results. For power law spectra, the recent WMAP results by themselves, favor a nearly scale invariant spectral index 
with n s — 0.99 ± 0.04, but when other large scale structure data is added slightly lower values of n s are favored 1 10]. Spergel et 
al 1 10] also explore more complicated, running spectral index models, for fitting the results from WMAP, other fine scale CMB 
experiments and large scale structure data. A recent study combining CMB and large scale structure data favors a scale invariant 
spectrum, with n s = 0.98 ± 0.02 with dn s /dln k = 0.003 ± O.Oli lllll . Slightly different set of parameters are derived by fl2il . 
when they combine the WMAP data with the SDSS results. 

One can relate the normalization constant A^ to the scalar normalization A used in CMBFAST and by the above authors. 
Verde et al |21] give A 2 (fc ) = (800tt 2 /T)^ 2 , where T = 2.725 x 10 6 /iK and k = 0.002 Mpc" 1 . For n = 1, and adopting 
a value A — 0.9, the best fit value for WMAP data alone, gives A$ ~ 3 x 10~ 5 , in agreement with the earlier COBE results. 
We can also relate A$ to the normalization of the matter power spectrum 8 2 H = k' i P(k)/(2n 2 ), evaluated at k = H . We have 
5h = (2/3)A ! / > (Di(a = 1) /SI™) 1 / 2 , where Di(a) is the growth factor fiUl . For a flat matter dominated model one would then 
get Sjj ~ 2 x 10 -5 consistent with earlier COBE results. 



B. The isocurvature mode 



The Sachs-Wolfe effect in theories which begin initially with isocurvature perturbations can be computed in an analogous 
manner. Suppose for example, one assumes that the universe has two dominant components, radiation with density pr and dark 
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matter with density p m . Then the total density perturbation will be St = {prSr + Pm8 m ) / (pr. + pm) — {$R + Z/<^m)/(l + V)- 
Here we have defined y = p m / Pr = aiv) / Q-iVeq) with Veq the epoch of matter-radiation equality. In such a model there is 
also an independent 'isocurvature' mode, where the initial curvature fluctuation is zero, but there are non-zero fluctuations in the 
'entropy per particle', S cx nn/n m cx T 3 /p m . This entropy fluctuation is characterized by a = 8S/S = 3Qo — 8 m . In terms of 
St and a, we have Sr — (<5t(1 + y) + j/<t) /(l + 3y/4). For such isocurvature initial conditions on super Hubble scales, the initial 
value of (T — Ci is preserved, and this drives the generation of curvature or potential perturbations in the radiation dominated era 
(see for example 1 20]). The resulting potential fluctuations freeze after matter domination, on super Hubble scales, and are given 
by = (Tj/5 with an associated density perturbation St = —20 = — 2c^/5. Also at matter domination, with y ^> 1, one has 
Sr -> (4/3)(5 T + cr). So 6 = Sr/4 -> (1/3) (-20 + 50) = and 8 + = 20. Then the associated CMBR anisotropies due 
to the monopole term in Eq. J37> is ai — (0o + Q)ji(kR*) — 2&ji(kR*). The associated dipole and quadrupole can again be 
neglected for super horizon scales. For the same amplitude of potential perturbations at the epoch of last scattering, isocurvature 
initial conditions therefore lead to 6 times larger temperature anisotropies on large scales. We have considered here only a two 
component system; for several components, several independent modes of perturbations can obtain with a variety of associated 
CMBR anisotropies 



C. The integrated Sachs-Wolfe effect 

If the potential were to change with time after decoupling, we see from Eq. J37i that further anisotropies can be generated at 
large angular scales. This effect is known as the integrated Sachs Wolfe effect, and typically arises in open universes, or in a flat 
universe with dark energy/cosmological constant, wherein the potential decays after the universe is dominated by curvature or 
dark energy respectively. The gravitational potential is also traced by other measures of large scale structure, There has therefore 
been considerable interest in checking whether there is a large angle cross-correlation between the temperature anisotropies 
(some of which will arise due the ISW effect) and other measures of large scale structure, with some tentative detections 1 23 ] . 

The Sachs-Wolfe and the ISW effects are dominant on scales larger than a few degrees, or I < 20 or so as schematically 
indicated in Figure 1 . In order to understand the smaller scale anisotropies we have to have study in greater detail the baryon- 
photon dynamics, to which we now turn. 



VII. THE BARYON-PHOTON DYNAMICS 



We have already derived the equation for the perturbed brightness for the photons. To complete the description of the baryon- 
photon system we have to also write down the continuity and Euler equations for the Baryons in the perturbed FRW universe. 
The continuity equation for the baryon density perturbation 8b is 
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(41) 



where term on the RHS takes account of the variation of the spatial volume due to the perturbed potential. In the baryon Euler 
equation, we include the force exerted by the radiation on the Baryons due to 7— e collisions. This force is most simply calculated 
as the negative of the rate of momentum density transfer to the photons by the electrons. The change in the momentum density 
of photons per unit conformal time to linear order in the perturbations is given by 
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(42) 



where (df/dr)) co u = c(f) is the change in photon momentum density due to collisions calculated in Section llVl The momentum 
transfer to the electrons will be negative of the value calculated in Eq. (142 \ . The radiative force density exerted on the electrons 
(and hence the baryons) by the radiation is then 
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(43) 



where F — J (dil/Air) i n is the first moment over photon directions of the fractional brightness. So the baryons feel a force 
due to the radiative flux F and a drag proportional to their velocity. The Euler equation for the baryons is then 



p B [v + Hv] = -p B V4> - Vp B + n e a T PRa 



(44) 



where Ti = (da/drf)/a and p B is the baryon pressure. The equations Eq. ( I26l l, Eq. J41> . Eq. (I44i together with an equation of 
state for the Baryon gas form the basic set of equations for the Baryon-Photon system. These equations can be solved to a good 
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approximation in the tight coupling limit, where we consider scales of the perturbations much larger than the photon mean-free 
path. This approximation is likely to be very accurate, before the recombination epoch, when matter is mostly in an ionized 
form. Note that the co-moving photon mean free path L 1 — (rieffra)" 1 grows to about Mpc scales just before the decoupling 
epoch. So the approximation kL 1 <C 1 (which corresponds to the limit I = kR* -C R* jL 1 ~ 10 4 ), should hold quite accurately 
for most scales of interests probed by CMBR anisotropics. One can then solve for the brightness perturbation iteratively. For 
this we first rewrite equation Eq. J26t as 



i = t H — h 4n • v - L n 



di 
dr) 



(45) 



(The repeated m index is assumed to be summed over). We can now write down the solution by iteration in powers of L 1 . We 
get 



■(0) ■ , l 2mY2m , , 

i = to H ^ ' 



(46) 



i« = j(o) _ jj 



drj 



Vz (0) + 8n • Vfi 



(47) 



Here and iW are iterative solutions to Eq. ( 145 1 giving i to the zeroth and first order in L 1 , respectively. (We will later 
consider iteration up to the second order when deriving Silk damping). 

Consider to begin with the effects of the baryon - photon tight coupling to the first nontrivial order given by Taking the 
zeroth moment of Eq. J47> . that is averaging both sides of Eq. (I47i over the all directions of the photon momentum, we get 



«o = «o 



where we have used the fact that J (dCl /4tt) 
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b, we then have 
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-V • v = 



(48) 



(49) 



This implies that the fractional perturbation to the photon number density Sna/nu = 39o satisfies the same equation as 8b- So 
initially adiabatic perturbations in the baryons, with 5g = 39q = (3/4)i5r initially, maintain this relation in the radiation era. 
The first moment, (that is multiplying Eq. ( I47l i by n and integrating over the directions of photon momenta) gives 
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The radiative force experienced by the baryons is then 
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So the Euler equation for the baryons, after substituting 6 = «o/4 + 4>, becomes 



[pB + ^Pr]v 



p B Hv = -[p B 



4 48 
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(50) 



(51) 



(52) 



We see therefore that in the tight coupling limit, the effect of Thomson scattering by radiation, to the leading order, is to add to 
the baryon Euler equation : (i) a radiation pressure gradient term with p ra d — Pi?40 o /3 = p^(5 7 /3, (ii) an extra inertia due to 
the radiation by adding a mass density (4/9^/3), to the inertial term in the LHS of Eq. d52t and to the gravitational force term 
on the RHS. When the radiation energy density and pressure dominate over that of matter, the baryon photon fluid, in the tight 
coupling limit, behaves as though its mass density is (pr + pr) and its pressure pn = pn/3 is due to radiation. The ratio of 
the inertia due to baryons and that due to radiation is given by R — Spb/^pr ~ 0.6(ilb/i 2 /0.02)(z/10 3 ) -1 . So baryon inertia 
cannot be neglected. (On the other hand the fluid pressure can be neglected compared to the radiation pressure, since the thermal 
speed in the fluid is much smaller than c/V3). 

On taking the time derivative of the continuity equation Eq. ( I49> . substituting for v from the Euler equation Eq. (I52> . and 
taking its Fourier transform, we get 
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We see that the baryon photon fluid can undergo acoustic oscillations, driven by the potential, and with an effective sound speed 
c s = 1/ \/3(I + R). If the baryon inertia were subdominant, with R — > 0, c s — > 1/ V3 which is the sound speed for a highly 
relativistic fluid. The baryon inertia leads to a reduction of c s from this extreme relativistic value. The oscillator equation 
Eq. \53\ can also be cast in a more suggestive form (cf.Eq. 16 of 0), 




We will use the solution of the oscillator equation Eq. (1531 to discuss the imprint of the acoustic waves on Ci . 



VIII. ACOUSTIC PEAKS 



The acoustic oscillations of the baryon-photon fluid lead to a rich structure of peaks and troughs in the CMB anisotropy power 
spectrum, on sub degree angular scales (or I > 100). To understand their basic features, let us look at an approximate solution 
of the oscillator equation Eq. J53I . To begin with let us neglect the slow variation of R with time, compared to the oscillation 
frequency kc s . Then we can rewrite Eq. ( I53l > as 

+ *»c?(6o + *) = + 20 (55) 

Also consider first modes which enter the Hubble radius in the matter dominated era, for which (d 2 §)/drj 2 w 0. Then the 
solution of Eq. d55i is 

6 + $ = A(k) cos kr s (i]) + B(k) sin kr s (rj) - i?$ (56) 

where r s (r)) = drf / a/3(1 + R) is called the 'sound horizon'. Note the sine and cosine oscillations will persist in the full 
solution but will have a slow damping due to a variable R. The — i?<I> term is the particular solution of the inhomogeneous 
equation. The effect of a non-zero R (called 'baryon loading') is to change the sound speed c s and also shift the zero point 
of the oscillations of the monopole (0q + $). One needs to specify initial conditions to fix A(k) and B(k) in Eq. d56l >. Note 
that as i] — > 0, for adiabatic or curvature perturbations, we already showed in Section IVT1 that 9o + $ — > 3>/3. This fixes 
A(k) — ($/3)(l + 3R). Also in the tight coupling limit, we have from Eq. ( I49K ikV = —dQ/drj. Using this relation, and 
noting from Section lvTl that for adiabatic initial conditions V — > as krj — > 0, fixes B(k) = 0. Imposing these initial conditions 
we have for modes which enter in the matter dominated era, 

$ 

©o + * = ~(l+3iZ) cosfcr s 
where we have again neglected the time variation of $ and R 
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(57) 



A. Radiation driving 



For modes which enter the Hubble radius during the radiation dominated era, one cannot neglect the variation in the grav- 
itational potential $. The comoving wavenumber k eq , corresponding to the Hubble radius at matter-radiation equality, is 
k eq — H(z eq )/(l + z eq ) — (2il m iJQZ eg ) 1 / 2 , and modes with k > k eq enter the Hubble radius in the radiation dominated 
era. During radiation domination, for the fluid which has an equation of state p = p/3, as would obtain for the tightly coupled 
baryon-photon fluid, the Einstein equations give $ + (4/77)$ + (fc 2 /3)$ = (dots as before denote derivatives with respect to 
conformal time) [20]. The solution for 'adiabatic' initial condition is 

3 

$>(k,rj) — - — [sinwfi — lot] cos ujrj] $j(fc) (58) 

where uo — k/^/3 = kc s is the frequency of the acoustic waves and $i(fc) the initial potential perturbation on super horizon 
scales. (Note that during radiation domination c s = l/v3). One sees that at early times on super horizon scales, with kr\ <C 1, 
$ — * where as once a mode enters the Hubble radius, the potential decays with time, going asymptotically as $>{k,rj) — > 
— (3 cos lot))$ i / (ujrf) 2 for krj 3> 1. This decay of $ causes extra driving of the acoustic oscillations for such modes. We can 
estimate the effect of this extra driving by directly solving for the associated density perturbation 5n — 5pji/pn using the 
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Einstein equations (cf. |H); 4nGa 2 Sp R = -fc 2 $ - - 3H 2 $, and 3H 2 = 8nGa 2 p R . For fcr; < 1 one gets S R -> -2$. ( , 
giving an initial value for the monopole (0 O + $) = S R /4 + $ — > 3>j/2. On the other hand, after a mode enters the Hubble 
radius, one has asymptotically, 5r — > — (2k 2 rf (k, rf) = 6$>i(k) cos(kc s r)) for 3> 1. So a mode which enters the Hubble 
radius early in the radiation dominated era has acoustic oscillations with 

e + $^e = <V4= !$i(fe)cos kc s n. (59) 

The amplitude of the oscillation, is therefore enhanced relative to a mode entering in the matter dominated era, by a factor 
(3$j/2)/(j|o/3j = 5, where we have used $o = (9/10)$j. This enhancement is referred to in the literature as 'radiation 
driving' l24l 12511 . The factor of 5 derived above gets modified to about 4, if we include the neutrino component 1 24] . It is 
also valid only in the asymptotic limit of very small scales and ignores the damping effect to be discussed below. Further the 
modes which are seen as the first few peaks in the C; spectrum have k/k eq larger than unity only by a modest factor, and so the 
enhancement is smaller. Nevertheless, the rise from the Sachs-Wolfe plateau of the C; versus I curve as I increases from a few 
10's to above 100 or so, as displayed in Figure 1, is dominated by this radiation driving effect. 

Note that due to the decay of the potential $, the baryon loading term i?$ in Eq. d57t is absent for modes which enter the 
Hubble radius well into radiation domination; so if one does see the effect of baryon loading in the C[s at higher /, this would 
be a firm evidence for the importance of a dark matter component in the universe (see below). 



B. Silk damping 



So far we have ignored the effects of departures from tight coupling. This departure introduces viscosity and heat conduction 
effects, and associated damping of the acoustic oscillations on small scales, worked out by Silk To calculate Silk damping 
effects, one needs to go the second order in L 1 . We give a detailed derivation, starting from the Boltzmann equation in Ap- 
pendix|B] In this derivation we neglect the anisotropy of the Thomson scattering, and also the effects of <f>. (The scales for which 
damping is important, enter in the radiation era, and so <\> decays as explained above). 

For plane wave solutions of the form, 



v = V exp(ifc • x 



Tdrj'); Go = Oo exp(ife • x 



Ydff) 



we derive the dispersion relation 
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(61) 



To the first order in L y , the baryon-photon acoustic waves suffer a damping, with the damping rate being larger for larger k or 
smaller wavelengths. This damping effect j^L is referred to in the literature as Silk damping, (If one takes into account the 
anisotropy of the Thomson scattering one gets 16/15 instead of 4/5 in the last factor above 1 27]). Silk damping introduces an 
exponential damping factor exp —(k/ko) 2 into the sine and cosine terms of Eq. d57i . where the damping scale kr> is determined 

by, 
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(62) 



(Also since modes with k > kjj, for which Silk damping becomes important, enter in the radiation dominated era and their 
potential $ has already decayed significantly; so the i?<l> term for such modes is not important). Since R grows to at most ~ 0.5 
by decoupling, the Silk damping scale fc^ 1 ~ [?7*L 7 (77*)] 1 / 2 by the last scattering epoch, or the geometric mean of the comoving 
photon mean free path and the the Hubble scale at last scattering. 



C. Putting it all together 

We can now put all the above ideas together to explicitly write Ci incorporating the baryon-photon oscillations. For scales 
much larger than the thickness of the LSS it suffices to use Eq. J37l > for a;, substituting the tight coupling expressions in Eq. j57l i. 
for 9q + <1> and \V . (The quadrupole term has negligible effect in the limit of tight coupling). The resulting a; is substituted into 
Eq. (|8} to compute Q. Then we have for the anisotropy power spectrum, 
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(63) 
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Note that Eq. d57i only describes accurately modes which enter in the matter dominated era. For modes which enter the Hubble 
radius during radiation domination, one has to take account of the k dependent enhancement due to radiation driving. Also for 
large k we have to take account of Silk damping. These effects can only be accurately incorporated in a numerical solution 
for C;. However many of the physical effects governing the properties of the acoustic peaks can be illustrated without such a 
detailed solution, keeping in mind that the co-efficients of the oscillatory terms will have an extra k dependence due to radiation 
driving and Silk damping. The fudge factor E(k) has been incorporated into Eq. (I63i to remind ourselves of the existence of 
these effects. 

It is also important to recall that ji(kR*) is a function sharply peaked at kR* ~ I. So for any given I, the k integral is 
dominated by modes which satisfy k ~ l/R*. On the other hand, the function j'Ax) is not as strongly peaked as ji(x) and has 
also a much smaller amplitude compared to ji ( see f° r example ll24l l28ll ). So the contribution from the doppler term (which 
contains j'JkR*)), is subdominant compared to the term depending on the temperature and potential (which contains ji(kR*)). 
We can now use Eq. (I63> to understand various features in the Ci spectrum. 

• The CMBR power spectrum, or C; has a series of peaks whenever the monopole term is maximum, that is when 
cos(fcr*) = ±1, where r* = r s (rj*) is the sound horizon at last scattering. This obtains for kr s {rf) = mr, where n 
is an integer; or for / ~ kR* = jiIa, where we define I a — ir(R* /r s (rj*)). These acoustic peaks were a clear theoretical 
prediction from early 70's 1 29, 30]; they used to be called doppler peaks, but note that the doppler term is subdominant 
compared to the temperature and potential contribution to Ci . The peak structure for a standard ACDM model is shown 
in Figure 1 . 

• The location of the first peak depends sensitively on the initial conditions, (isocurvature or adiabatic) and also most 
importantly on the curvature of the universe. The current observations favor a fiat universe. For a fiat geometry, the 
location of the first peak can be used to measure the age of the universe. 

• For isocurvature initial conditions the monopole term would have sin(fcr*), which would be maximum at kr* = (2m + 
1)tt/2, where m = 0, 1, 2... The peak at kr* — n/2 is generally hidden. The first prominent peak for isocurvature initial 
condition is at kr* = 3tt/2, and so occurs at larger I than for adiabatic perturbations; present observations favor adiabatic 
initial conditions. 

• Almost independent of the initial conditions the spacing between the peaks is ~ I a. 

• Due to a non zero baryon density, that is a non zero R, the peaks are larger when cos(kr* s ) is negative, since in this case, 
the cosine term and the — 3i? term in Eq. d63i add. Due to this effect of 'baryon loading', the odd peaks, with n = 1,3,.. 
have larger amplitudes than the even peaks with n = 2,4, ... 

• The radiation driving effect, as we explained earlier causes the C; curve to rise above the Sach- Wolfe signal for I values 
corresponding to the acoustic peaks (cf. Figure 1). Note also that the i?$ term would be absent, if the the scale cor- 
responding to a given peak enters during radiation domination, such that the potential $ has decayed by the epoch rj*. 
Indeed the observed existence of a 3'rd peak almost comparable in height to the 2'nd is an indication of the importance of 
(dark) matter in the universe. 

• Silk damping cuts of the C; spectrum exponentially beyond I ~ kr>R* ~ 1500 (cf. Figure 1). There is also damping 
of the Cf spectrum due to the finite thickness of the last scattering surface. The scales for both damping are similar. The 
decline in C; due to both effects has been parametrized by a cxp[— (I/Id)" 10 ] factor, where Id = kpR*, and tud ~ 1.2 

We also mention some of the other consequences of the varying gravitational potential, for the C; spectrum. 

• Th effects of a varying gravitational potential lead to the ISW effect as mentioned earlier. This can operate both after 
last scattering and during the period of recombination. In an universe which is at present dominated by dark energy, the 
potential associated with sub horizon scales decay after dark energy domination. The resulting increase in C\ leads to the 
upturn for I < 10, from the Sachs-Wolfe plateau seen in Figure 1. 

• There is also an early ISW effect for modes which enter the Hubble radius in the radiation dominated era. However due to 
the e~ T factor multiplying $ in Eq. A36I . this contributes to C; only for those modes whose potential's decay just before 
last scattering. The early ISW causes an increase in C; for such modes. 

• The early ISW effect partly fills in the rise to the first peak and leads to a shift in the location of the first acoustic peak to 
a lower I < I a ltl5ll . Also for modes with k > k eq , entering the Hubble radius in the radiation era, the decaying potential 
leads to a difference between the exact solution to Eq. J53i from the approximate solution given by Eq. d56i . This leads to 
a further shift in the location of the acoustic peaks, to lower I iflSl. l31fl . Finally, jf(x) has a peak at slightly smaller I than 
I = x. All these effects lead to a shift of the peak location to an I value lower than I — uIa, by ~ 25% or so, which can 
only be calibrated by numerical solution 113 ill (see below). 
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Note that we can use both the location and the relative heights of the acoustic peaks as a sensitive probe of the cosmological 
parameters, an issue to which we now turn. 



D. The acoustic peaks and cosmological parameters 



The cosmological parameters which have been constrained include the curvature of the universe or the total energy density 
fix, the baryon density io h = fl^h 2 , dark matter density u> m = £l m h 2 (which is predominantly thought to be cold dark matter), 
and the slope of the primordial power spectrum n s . We outline some of these ideas, following mainly Hu et al 1I3TI1 and the post 
WMAP analysis of Page et q/IOl. 



1. The location of the acoustic peaks 



For the fiat matter dominated universe, the conformal time rj oc a 1 / 2 oc (1 + z)~ 12 . If we neglect the effect of baryons, 
c s = 1/V3 and r* = rj*/\/3. Also R* = t] Q — rj*, and so the acoustic scale Ia = V3n(r)o — V*)/v* ~ V3^(vo/v*) = 
172(z*/10 3 ) 1 / 2 . We therefore expect the first acoustic peak around this value. It is however important to also take account 
of the radiation and baryon densities before decoupling. Radiation density increases the expansion rate and the baryon density 
decreases the sound speed and so r* gets altered (cf. Eq 2 in ll33ln 
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(64) 



Here r* = PR{tf) / p m {rf) ~ 0.3(w m /0.14)~ 1 (z*/10 3 ). Also for a universe with non-zero curvature, in determining the 
mapping between I and k, its necessary to replace the comoving angular diameter distance d,A = R* = T)o — V* corresponding 
to a flat universe, by cLa applicable to a general cosmology. This is given by, Jl7ll24ll 
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(65) 



Here Q,y is the ratio of the dark energy to the critical energy density, and w the dark energy equation of state parameter (w = — 1 
for the cosmological constant). For a flat ACDM cosmology with = 0.73 and £It = 1 one gets d ~ 0.89. Using 
I a = ndA/r*, we see that the u> m dependence cancels out and 
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(66) 



Note that for a flat universe ignoring the effect of baryons and radiation, one then gets I a ~ 172, as before. But with ujb = 0.02, 
u> m = 0.14, even for a flat universe, I a ~ 300 and so is much larger. We note from Eq. i66\ that the acoustic scale is most 
sensitive to the value of £It, the total density parameter. 

Further, the location of the first peak is shifted from I a because of the effects of potential decay (as described above), which 
becomes important for modes with k >k &q , entering the Hubble radius during radiation era. The comoving wavenumber k, 
corresponds to I = l eq = k eq dA, where [24] 
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(67) 



One needs to work out the exact shift numerically; For a scale invariant model, with n s = 1 and to h = 0.02, Hu et al 1 3 JJ] give 
In = Ia(it> — VO where ip ~ 0.267(r*/0.3) 01 , and for better accuracy one replaces 0.267 with 0.24 for I2 and 0.35 for ^3. For 
example, for a flat ACDM cosmology with u) h = 0.02, u m = 0.14, 57t = 1, w = — 1, and taking account of the phase shift, the 
first peak is predicted to be located at 1% w 220. For the WMAP data, the measured value of the l\ — 220.1 ± 0.8. So the data is 
indeed consistent with such a flat universe. (The peaks also get affected mildly by the tilt in the power spectrum from n s = 1). 

However one should caution that l\ alone does not determine the geometry; one needs some idea of fi TO and flf, which can 
be got from the full WMAP data. There still remains potential degeneracies 1421 l43l 14411 . whereby the peak location can be left 
unchanged by simultaneous variation in tt m — h space and f2 m — Qji space. If one imposes h > 0.5 as seems very reasonable, 
one infers 0.98 < < 1.08 (95% confidence level) 1 10]. For the HST Key project measurement of Hq as a prior, one gets 
D.T = 1-02 ± 0.02. The observations strongly favor a flat universe. Also from the inferred values of ust and u> m from the full 
data, one gets an acoustic scale Ia ~ 300. If one assumes a flat universe, it turns out that that the position of the first peak is 
directly correlated with the age of the universe. The WMAP data gives t = 13.6 ± 2 vr for the ACDM model Qj| 
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Finally, the whole Ci spectrum is damped strongly beyond the scale Ijj = kjjR* . Numerically, we have from Hu et al 

2240 d ( z* 

D ~ [(i + r *)i/2_( r *)i/2]i/2 V103 

For the ACDM model with WMAP parameters, one gets Id ~ 1470. The damping scale shows a much stronger dependence on 
u>b and the redshift z* compared to The small angular scale experiments like the Cosmic Background Imager (CBI) |36] do 
find evidence for such a damping. 
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(68) 



2. Peak heights 

The heights of the different peaks, can also be used to infer cosmological parameters. We define the height of the first peak as 
HH, Hi = {AT h I AT 10 ) 2 , that giving its amplitude relative to the power at I = 10. (For the WMAP data the height of the first 
peak is AT^ = 74.7 ± 0.5/xK). Hi increases if (a) uj m decreases (because radiation driving is more effective a lower matter 
density), (b) if u>b increases (due to the baryon loading) (c) if one has a lower Q,& or higher f2y (because then the integrated 
Sachs Wolfe effect is smaller which decreases ATio). Further Hi can decrease if one has re-ionization (since a fraction r of 
photons are re-scattered), or if one has a contribution from tensor fluctuations (tensors will contribute to Sachs-Wolfe effect but 
not to acoustic oscillations). Since Hi depends on several effects, there is no simple fitting formula; around ACDM 113 ill have 
given a crude formula for its variation with various parameters. 

The height of the second peak is defined relative to the first, as H2 = (AT; 2 /AT; 1 ) 2 . This ratio is insensitive to reionization 
or to the overall amplitude of the power spectrum since these scale both peaks by the same amount. The dependence on oj m is 
also weak because radiation driving roughly affects both peaks similarly. Hi is most sensitive to the baryon density uif,, since 
baryon loading increases the first peak relative to the second. It is also sensitive to any tilt in the spectrum, away from n s = 1. 
From fitting to a grid of spectra using CMBFAST 1 19], one has 1 33] 

H 2 = 0.0264w b ~°- 762 (2.42)" s - 1 x exp[-0.476 ln(25.5w b + 1.84w m ) 2 ] (69) 



= 0.88An s - 0.67^ + 0.039^ (70) 

H2 LUb W m 

For the WMAP data, H2 = 0.426 ± 0.015. For a fixed uj m the first two terms of Eq. d70t quantifies the degeneracy in the ujb — n s 
plane. 

The height of the third peak increases as u)b increases (baryon loading). The ratio H3 = (A7] 3 / AX^) 2 is most sensitive to 
n s or any departures from scale invariance, because of the long / baseline. Hu et al ll3lll give 

H= 2.17^ 59 (3.6)"-- 1 

3 [1 + (w fc /0.044) 2 ][l + 1.63(1 - w b /0.071)w m ] 



^ = 1.28An s - 0.39^ + 0.46^ (72) 

Hz ujb u> m 

These dependencies are accurate to few percent levels for variation around the WMAP inferred parameters l33ll . WMAP does 
not yet clearly measure the third peak, but from previous compilations |45], Page et al estimate Hj, = 0.42 ± 0.08. Note that if 
n s is fixed, u>b is well constrained by H2 and then u> m from H3. For more details we refer the reader to I33tl . We show in Figure 
2 a set of C/ versus I curves, generated using CMBFAST, which illustrate the sensitivity of the CMBR to the cosmological 
parameters discussed above. 



E. Other sources of CMB anisotropics 



So far we have concentrated on the primary temperature anisotropics generated at the LSS; a number of effects can generate 
additional anisotropies after recombination, generally referred to as 'secondary anisotropics'. We do not discuss these in detail; 
for an extensive review see 111 411 . Of the gravitational secondaries, we have already discussed the ISW effect arising from 
the changing gravitational potential. This effect is also important if there are tensor metric perturbations, say due to stochastic 
gravitational waves generated during inflation 13411 - Another important gravitational secondary arises due to gravitational lensing 
(cf. 11411 and references therein). Scattering effects due to free electrons along the line of sight can also produce a number of 
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effects. The electrons can arise in collapsed objects like clusters or due to re-ionization of the universe. We discuss the effects 
of re-ionization later below. The scattering of the CMB due to the ionized electrons in clusters of galaxies was first discussed by 
Sunyaev and Zeldovich (SZ) 1 35 ] . The SZ effect generates power below the damping tail in the C; spectrum, at a level which 
depends on the normalization of the density power spectrum, erg. (Here erg is the RMS density contrast when the density field is 
smoothed over a 'top hat' sphere of radius 8/i -1 Mpc). Recently a significant excess power was detected by the CBI experiment 
136], at small angular scales ( / > 2000) at a level of ~ 355(^ K) 2 . This can arise from the SZ effect, but requires a somewhat 
large er 8 ~ 1 (cf. l37ln . larger than values previously assumed. Alternatively the CBI result may point to new physics; it has been 
suggested for example that primordial magnetic fields can be a significant contributor to the power at large I | 38, 39]. Primordial 
tangled magnetic fields generate vortical ( Alfven wave mode) perturbations which lead to temperature anisotropics due to the 
doppler effect. They also survive Silk damping on much smaller scales compared to the scalar modes 1 40, 41]. The test of 
whether the CBI excess is indeed produced by the SZ effect, will come from the spectral dependence of the excess power (if it 
is due to the SZ effect, there should be such a spectral dependence), and measurements of polarization on these small angular 
scales (see below). There are several other interesting gravitational and scattering secondaries which can generate temperature 
anisotropics, and we refer the interested reader to the excellent review 0. 

IX. POLARIZATION OF THE CMBR 

A. The origin of CMB polarization 

It was realized quite soon after the discovery of the CMB that it can get polarized |46]. Polarization of the CMBR arises due 
to Thomson scattering of the photons and the electrons, basically because the Thomson cross section is polarization dependent. 
We used in earlier sections the cross section relevant for unpolarized light, ignoring the small effects of polarization on the 
temperature evolution. Scattering of radiation which is isotropic or even one which has a dipole asymmetry is however not 
capable of producing polarization. The incoming radiation needs to have a quadrupole anisotropy. The general features of 
CMBR polarization are discussed in detail in some excellent reviews |47, 48]. Note that in the tight coupling limit, the radiation 
field is isotropic in the fluid rest frame, and can have at most a dipole anisotropy in the frame in which the fluid moves. The 
quadrupole anisotropy is zero. However to the next order, departures from tight coupling, due to a finite photon mean free path, 
in the presence of velocity gradients, can generate a small quadrupole anisotropy. 

A qualitative argument is as follows 1 16]: The last scattering electron (say O at x$) sees radiation from the 'last but one 
scattering' electron (P), roughly a photon mean free path (L 7 ) away, say at a location x = xq + L^n. Here n is the direction 
from O to P. The velocity of the baryon-photon fluid at P is Vi(x) as Uj(xo) + L- l rijdjVi(xo). Due to the Doppler shift, the 
temperature seen by O is ST(xq, n)/T ~ rii[vi(x) — Uj(cco)] = L 1 riinjdjVi(xo). This is quadratic in n and so corresponds to 
a quadrupole anisotropy as seen by the last scattering electron. The Thomson scattering of this quadrupole anisotropy can lead 
to polarization of the CMBR. The fractional polarization anisotropy generated is Ap ~ kL y V. 

One complication is that L 1 grows rapidly as photons and baryons decouple during recombination. An approximate estimate 
of its effect, would be to weigh the polarization amplitude derived above, with the probability of last scattering at a given epoch 
described by the visibility function. Note that the visibility function goes as re~ T , where f = 1/L 7 . So during the tight coupling 
evolution, the L 1 factor cancels out and after the weighting one gets instead Ap ~ kSrj*V, where 8rj* is the width of the LSS. 
So the effective photon mean free path generating quadrupole anisotropy and hence polarization of the CMB becomes 6r/*, the 
average distance photons travel between their last and last but one scattering, during decoupling. Such an estimate is verified in 
a more careful calculation i49Tl . 

B. Describing CMBR polarization 

There is another complication that has to be handled when dealing with polarization, the fact that polarization is not a scalar 
quantity. It is conventional to describe polarization in terms of the Stokes parameters, /, Q, U and V, where / is the total 
intensity, whose perturbed version was called i above, and discussed extensively in earlier sections. For a quasi-monochromatic 
wave, propagating in the z-direction, we can describe the electric field at any point in space as E x — a x (t) cosfwoi — d x (t)] 
and E y = a y (t) cos[uiot — y (t)], where the amplitudes a x , a y and the phases 9 X , 8 y vary slowly in time, compared to u>q . 
The stokes parameters are defined as the time averages: / = (a 2 ,) + (a 2 ), Q = (a 2 ) — (a 2 ), U — (2a x a y cos(9 x — y )), 

V = (2a x a y sm(0 x — 6 y )). Unpolarized light has Q = U = V = 0. The parameters Q and U describe linear polarization, while 

V describes circular polarization. At the zeroth order the CMB is unpolarized and its small polarization is expected to arise as 
explained above due to Thomson scattering. This does not produce circular polarization and so one can set V = 0. 

Note that under a rotation of the x and y axis through an angle ip, the parameters / and V are invariant but (Q ± iU) 1 = 
e T2iif> _|_ jjj-j So QzLiU transform as a spin 2 Tensor under rotation of the co-ordinate axis. The standard spherical harmonics 
do not provide the appropriate basis for its Fourier expansion on the sky. One then adopts the following approach to this problem 
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150115111 : construct scalars under rotation from Q ± iU by using spin-lowering (d~) and spin-raising (d + ) operators, and then 
make a standard Yi m expansion. Or alternatively construct tensor ('spin' weighted) spherical harmonics, ±2^m by operating on 
the Yi m 's twice with spin-raising or lowering operators, and then expand 



a±2,lm (±2^™) 



(73) 



in this basis. Alternatively a±2j m can also be thought of as the Y\ m expansion co-efficients of the spin zero quantities, (d~ ) 2 (Q+ 
iU) and (d + ) 2 (Q — iU) respectively, apart from an / dependent normalization factor. The explicit expressions for the raising 
and lowering operators, the spin weighted harmonics, and the expansions in terms of these are given in |50S]. For example we 
can write, 



a±2,, 



±2 



Y l * m (n))(Q±iU)(n) 





~"7 


l(l-2)\_ 





<mY i ; n (n)[(dTf(Q±iU)(n)] 



(74) 



Since a±2 t i m are expansion co-efficients of scalar quantities under rotation, they can be used to characterize the polarization 
on the sky in an 'invariant' manner. More convenient is to use the linear combinations, a,E,im = — (»2,Zm + »-2,2m)/2 and 
a B,im — K a 2,im — a -2,2m)/2 HO EH- and the associated real space polarization fields; E(n) = YLim a E,imYi m (n) and 
B{n) — J2im a B,imYi m (n). The E and B fields specify the polarization field (Q and U) completely, are invariant under 
rotation (just like the temperature 8(n)) and have definite parity. Under a parity transformation, E remains invariant while B 
changes sign 1 52] . The convenience of the E-B split comes from the fact that scalar perturbations do not produce any B type 
polarization. An alternative way of thinking about the E and B split is that they are the gradient and curl type components of 
the polarization tensor |51]. More details of these fascinating but somewhat complicated ideas can be got from the two seminal 
papers on the subject 1 50. Dill . 

In order to describe the statistics of CMBR anisotropies fully, including its polarization, we have to now consider not only C; 
due to the temperature anisotropy 0, but also corresponding power spectra of E, B and the cross correlation between 9 and E. 
Note that the cross correlation between B and 9, and B and E vanish if there are no parity violating effects. Since E and B 
are rotationally invariant quantities, we can define the power spectra Cf, Cf and Cf E in an analogous way to the temperature 
power spectrum. We now turn to their computation. 



C. Computing the polarization power spectrum 

We focus on scalar perturbations. In this case for any given Fourier mode k, one can define a co-ordinate system with 
k || z, and for each plane wave, treat the Thomson scattering as the radiative transport through a plane parallel medium. It 
turns out that only Stokes Q is generated in this frame because of azimuthal symmetry, and its amplitude depends only on 
(j, = n ■ k. The stokes parameter U = in this frame, for each k mode. Because U = and Q is only a function of /i, one has 
(d~) 2 (Q + iU) = (d + ) 2 (Q — iU). (From the explicit form of the spin-raising/lowering operators give in 1 50], it can be checked 
that (9~) 2 (/(/i)) = (<9 + ) 2 (/(/i)) for any azimuthally symmetric function which depends only on /i. Second since U = 0, we 
have (Q + iU) = (Q — it/)). Therefore 0,2. / m = a-2.im, and we have a,B,im — for such scalar perturbations. 

The Boltzmann equation including polarization is given by a number of authors (see for example 1321.15 51 56]). We will simply 
quote the result, got using the detailed treatment by 1331 . We have 

Cf = lj *fc 3 (|af(Mo)| 2 ); af{k, m )=-J d V g( V0 ,r,)^ 

where we have expressed the quadrupole source for the polarization anisotropy, P = (O2 — \/6i?2)/10, by its tight coupling 
limit P = 2kL 1 V/9 (see Bi t. As argued on qualitative grounds above, polarization is sourced by the velocity differences of 
the fluid, over a photon mean free path (i.e. kVL 7 ). Once again the spherical Bessel function ji (kArj) in Eq. d75i will at a given 
I, pick out (on k integration) a scale k ~ 1/ At] at around last scattering, while the visibility function g weighs the contribution 
at any time 77 by the probability of last scattering from that epoch. 

Suppose we wish to estimate the polarization anisotropy on physical scales much bigger than the thickness of the last scattering 
surface, or I ~ kR* < 1000 or so. As we explained earlier, the visibility function goes as fe~ T whereas the polarization source 
is kL^V/3 = (fcy/3)(f) _1 , and so in their product, f cancels and only e~ T (kV/3) would survive. The integral over r\ in 
Eq. ( 1751 . would be nonzero only for a range of epochs of order the width Srf of the LSS. (Note that just after recombination, 
the tight coupling expression cannot be used; however there is also no polarization for 77 > 77* because there is negligible further 
Thomson scattering). So one expects a contribution of order kVSij* /3 in doing this integral, apart from an evaluation of the 
other terms at 77*. A more rigorous analysis, following the time evolution of the polarization source term, gives a further factor 



(1 + 2)1 



a-2)\ 



1 MkArj) 
(kArj) 2 



(75) 
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of ~ 1/2 reduction, if 8rf, is defined as the Gaussian width of the visibility function i49ll . Making such an approximation, and 
putting in the tight coupling expression for the velocity of the photon-baryon fluid, gives 

of (Mo) = -^f- *o E(k)(l + 3R)c a smkri* s U^j 3l (kR*) (76) 

Note that again the k integral to find Cf will pick out values of kR* ~ I. We can infer a number of features of the polarization 
from the above: 

• The magnitude of the polarization anisotropy, is of order AP ~ 0.6(fc<5r;*)(<I>o/3) = 0.6 l(5r]* /R*)($o/3), where we 
have taken R ~ 0.6. Adopting Sr/* ~ 10h~ 1 Mpc and R* ~ 10 4 /i _1 Mpc, we get at I ~ 100, a polarization anisotropy 
about 6% of the Sachs-Wolfe contribution (or about 2/iK). The amplitude rises with I, but at large I > Id the Silk 
damping cuts off the baryon-photon velocity, and so the polarization gets cut off as e~^^ ln ^ m , say. This gives a maximum 
contribution at I < Id depending on the m, with peak amplitude of order 10% of the peak temperature anisotropics. These 
order of magnitude estimates are borne out by the more detailed numerical integration using CMBFAST shown in Figure 
1. 

• The acoustic oscillations of the baryon-photon fluid velocity imprints such oscillations also on the polarization. The 
polarization will display peaks when sin(fcr*) = ±1, or for fcr* = (2ri + l)(7r/2), with n = 0, 1, .., corresponding to 
I ~ (2n + 1)?a/2. These peaks are out of phase with the temperature acoustic peaks, as they arise due to the velocity, and 
they are sharper (since for temperature there is a partial filling in of the troughs by the velocity contribution). 

• Both the polarization and the temperature depend on the potential $, and so one expects a significant cross correlation 
power Cj E . Further, the ji term does not significantly correlate with j', term in the fc-integral for Cj E . So this cross 
correlation will be dominated by the product of the temperature monopole with a cos(fcr* ) dependence and the polarization 
(of E type) with a sin(fcr*) dependence. The peaks of Cj E will then occur when sin(fcr*) cos(fcr*) oc sin(2fcr*) = ±1, 
or when kr* = (2n + l)(jr/4), with n = 0, 1, .., corresponding to I ~ (2n + l)/^/4. So Cj E has oscillations at twice 
the frequency compared to the temperature or polarization. There will be shifts in the exact location of the Cf and Cf E 
peaks, as for the temperature. 

• The E type polarization has been detected at a ha level by the Degree Angular scale Interferometer (DASI) at I values 
of a few hundred, and at a level consistent with the expectations from the detected temperature anisotropy 1 53 ] . The CBI 
experiment has also detected the E type polarization, with the peaks in the polarization spectrum showing the expected 
phase shifts compared to the peaks of the temperature spectrum 15411 - The TE cross correlation was detected at 95% 
significance by DASI, but there is no evidence of any B type polarization. The cross correlation has also been detected by 
WMAP. The WMAP experiment has released results on C, , although not on Cf. WMAP detects significant negative 
Cf E , at I ~ 150 and a positive 'peak' at I ~ 300. The existence of such an anti-correlation between temperature and 
polarization is an indication that there exist 'super-Hubble' scale fluctuations on the LSS. This is interpreted as strong 
evidence for inflation type models, since models which involve seeds (like cosmic strings) can produce super Hubble scale 
temperature fluctuations (due the ISW type effects) but not the observed anti-correlation in Cf E . 



D. B type polarization 



So far we have emphasized the E type polarization, as scalar modes do not produce the B type signal. However models of 
inflation which are thought to generate the scalar perturbations, can also generate a stochastic background of gravitational waves. 
These tensor perturbations and the CMBR anisotropy that they generate has also been studied in detail |57], although we will not 
do so here. Their effects are best separated from the scalar mode signals, by the fact that Tensors also lead to B type polarization 
anisotropy 11501. Bill . The temperature contribution from tensors is flat roughly upto I ~ 100 after which it rapidly falls off. The 
polarization contribution, produced at recombination, peaks at I ~ 100. The peak amplitude of the signal is however expected 
to be quite small in general with (1(1 + l)Cf /27r) 1/2 — 0.1/j,K(E in f /2 x 10 16 GeV) 2 , where E mf is the energy scale of 
inflation. II 611. (An Ei n f ~ 2 x 10 16 GeV corresponds to the ratio of the I = 2 contribution due to tensors compared to scalar, 
T/S « 0.1). One of the prime motivations for measuring polarization with great sensitivity is to try and detect the contribution 
from stochastic gravitational waves. The B type anisotropy can also arise due to gravitational lensing of the CMB, even if one 
had only E type polarization arising from the recombination epoch 1 58 ] . This could set the ultimate limitation for detecting the 
B mode from gravity waves. Another interesting source for B type polarization are vector modes, arising perhaps due to tangled 
magnetic fields generated in the early universe I^l5l*l60ll6lll . or even present in the initial conditions |62]. Indeed if there 
were helical primordial magnetic fields, at the LSS, parity invariance can be broken and one can even generate T — B cross 
correlations |63||. 
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E. Reionization and CMBR polarization 

One of the surprises in the WMAP results was the detection of a significant excess cross-correlation power Cf E at low 
I, over and above that expected if polarization was only generated at recombination. This can be interpreted as due to the 
effects of re-ionization, But one seems to need a significantly higher optical depth to the re-ionized electrons r ri ~ 0.17, and 
a correspondingly high redshift for reionization z r i ~ 17. The probes and models of the high redshift intergalactic medium, 
including the use of the CMBR as a probe of re-ionization is discussed more fully elsewhere in this volume by Sethi 16411. We 
make a few qualitative remarks. 

First, note that if photons are re scattered, due to electrons produced in re-ionization, the visibility function will have 2 peaks; 
one narrow peak around recombination, and a broader peak around the re-ionization epoch (cf. Figure 2 in 1 64]), which depends 
on the exact re-ionization history. The probability for last scattering around the usual LSS will diminish by a multiplicative 
factor e~ r ", where r r i is the optical depth for electron scattering to the re-ionization epoch. At the same time new temperature 
and polarization anisotropies get generated. The most important effect is that Thomson scattering by electrons generated during 
re-ionization, produces additional polarization. Note that the quadrupole anisotropy at the the re-ionization epoch is likely to 
be much larger than at recombination, simply because the monopole can free stream to generate a significant quadrupole at 
the new LSS. At re-ionization redshifts close enough to the observer, the relevant monopole becomes the Sachs-Wolfe value 
Go + $ = $o/3- The quadrupole at the re-ionized epoch r) r i can then be simply estimated by replacing 770 in Eq. (I37> . by 
rjri. One gets 62(^,77,.^) = a2(k,i] r i) — ($o/3)j2(fc(»7rt — V*))- Note that this does not have the kL 1 suppression factor, 
which obtains around recombination. Also E% in the polarization source term P above is negligible. In evaluating the E type 
polarization arising from the re-ionization, one can substitute the resulting P — 02/10 in Eq. |75i instead of P = 2kL 7 V/9; 
for the range of r\ where scattering by electrons generated due to re-ionization is important. 

The resulting re-ionization contribution can be best calculated numerically, for example using CMBFAST, as illustrated by 
Sethi (this volume). But the scale where the peak in the power spectra can be estimated noting that af will involve the product 
3i{k{ri r i — T]*))ji(k(r]Q — r) ri )), which contributes to the fc-integral dominantly when both k(r\ ri — 77*) ~ 2 and k(ri — rj ri ) ~ I. 
This implies that the re-ionization contribution to E type polarization peaks at I ~ 2 (770 — r] r i)/ (j] r i — 77* ) ~ 10 for the parameters 
appropriate for a ACDM cosmology and a z ri ~ 20. This scale basically reflects the angle subtended by the Hubble radius at 
re-ionization. One has to also take account of the damping due to the large width of the LSS at re-ionization, which will shift 
the peak to smaller I. 

The k integral which determines Cf E , involves the product j2(k(f] r i — r)*))ji(k(r)o — f] r i))ji(kR*), the last ji(kR*) coming 
from the temperature contribution from the usual LSS. Note that 770 is much bigger than both r\ r i and 77*, and the two ji factors 
will re-inforce each other for small I. The cross correlation peak will occur at an I similar to the peak in Cf . The indication from 
the WMAP data for significant optical depth from re-ionization is not easy to explain (cf. 1 65 ]). If the preliminary WMAP result 
continues to firm up with subsequent years data, it will set very strong constraints on the star and active galaxy formation at high 
redshift. It may be also worth exploring new physical alternatives. For example Ref ll66ll . explores the possibility that tangled 
magnetic fields generated in the early universe could form subgalatic objects at high enough redshifts to impact significantly on 
re-ionization. Note that if there is significant optical depth to re-ionization, then inhomogeneities at the new LSS can lead to 
new secondary sources of both temperature 1 67] and polarization anisotropies |68|. Eventually the detailed measurement of the 
polarization signals, could be a very effective probe of the reionization history of the universe ll69il . 

X. CONCLUDING REMARKS 

In this review we have tried to emphasize the physics behind the generation of CMBR anisotropies. We have given details 
of the computation of the primary temperature anisotropies, and also indicated the relevant issues for polarization. Our aim is 
more to introduce the budding cosmologist to the well known (and reviewed) techniques used to calculate the CMB anisotropies, 
rather than provide an extensive survey of observations and results. Of course, it is the existence of very good observational 
data that makes the effort worthwhile. Clearly the CMB is and will continue to be a major tool to probe structure formation and 
cosmology. We have already learned a great deal from the detailed observations of the degree and sub degree scale temperature 
anisotropies, particularly the acoustic peaks. The exploration of small angular scale anisotropies is just at a beginning stage 
and holds the promise of revealing a wealth of information, on the gastrophysics of structure formation. The future lies in also 
studying in detail the polarization of the CMBR. Already WMAP results have revealed a surprisingly large redshift for the re- 
ionization of the universe. Polarization will also be a crucial probe of the presence of gravitational waves. We can expect in the 
years to come much more information on cosmology from WMAP, future missions like PLANK and other CMB experiments, 
with the possibility of more surprises! 
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APPENDIX A: THE COLLISION TERM: DETAILS 

Consider the integral over the collision term on the RHS of Eq. ( 1251 . We have J p 3 dpc(f) — n e <jT(Ai + A 2 ) where 

A i = J P 3 dp^[f(P, n') - f(p, n)]; A 2 = J p 3 dp~P 2 (h ■ h')[f(p, n') - f(p, n)] (Al) 

are respectively, the isotropic and anisotropic contribution to the collision term. From the invariance of the scalar piU 1 , where u l 
is the four velocity corresponding to the bulk motion of the electron (Baryonic) fluid we can show that 



p = a(l + n.v)(l — <j))p 



(A2) 



(We have used here the fact that in the fluid rest frame the components of u l — (1, 0, 0, 0), while u l = (j v / y/goo, Jvv/ \/\g^jf\) 
with 7 t , = (1 - v 2 /c 2 )~ 1/2 )- We split A X = I X - I 2 with 



h= j p 3 dp^-f{p,n>) I, 



3 , dSV j._ _.. 
P dp—f(p, n 
An 



(A3) 



For evaluating I 2 we use the fact that / is a scalar, that is f(p, n) = f(p, n). Also the integrand of I 2 does not depend on n'. 
So we have 1% = J p 3 dpf(p, n) = (p^a 4 /47r)[l + i]. For evaluating I\ we stay in the initial electron rest frame and transform 
the integral overp to one overp using Eq. &A2\ . We get I\ — (a 4 /47r)(l — 4(/))(l + 4n ■ w)p where we have used the fact that 
p = J d 3 ppf(p) is the energy density of radiation in the fluid rest frame. Using the invariance of Tlutu k , and from the fact that 
the components of both u l and which involve one spatial index are of order v/c, we can check that p = p + 0(v 2 /c 2 ). Since, 
p = p R {l + 40)(1 + i ), to linear order l x = (a 4 p R /ATi)(l + 4n • v + i ). So 



Ai = —Pr [io + An-v- i] 
To simplify A 2 , we use the addition theorem for spherical harmonics to write 
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m=-2 



Y 2m (n) 
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p 3 dp dfi' Y* m (n')[m n') ~ f(P, «)] = 



a PR 1 
4tt 10 



(A4) 



(A5) 



where i 2m = J dVlY 2m i(x 1 rj, n). In evaluating A 2 we have used the fact that the term f(p, n) does not contribute to the integral 
over dVif. Also writing f(p, n') = f(p, n') = fb(p) + fi(p, n'), the /& term gives zero contribution. And since f x is already 
first order in perturbations, we can evaluate A 2 by replacing n and n' by their unbarred values (these will differ only by terms 
of order v/c and the difference when multiplied by i 2m will not contribute to the first order). 

Finally, we also need to evaluate dr/drj. Since A\ and A 2 are already of first order, we need to evaluate this term only to 
zero'th order, to write down the equation for the perturbed brightness. We have dr jdr\ = u°/u° = a to the leading order. The 
perturbed brightness equation (1261 . given in the main text, is got from Eq. d25l >. Eq. iAAi and Eq. JA5l l. 



APPENDIX B: SILK DAMPING: DETAILS 

We have given in the main text the iterative solution to Eq. ( 145 \ to the first order in L 1 . To derive Silk damping one needs to 
go to the second order iteration, 
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As mentioned in the text, we neglect the anisotropy of the Thomson scattering, and also the effects of the gravitational potential 
<f). Taking the zeroth moment of Eq. dBll , we get 
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So to the next order in L 7 , Eq. J49> is modified to 
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Similarly, taking the first moment of Eq. (IB 1 i the Euler equation Eq. ( 1521 gets modified to 
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Here we have used the relation J (e2Q/47r)njn,-7i/.ri; = + ^ife^j; + 6ij6jy]/15. (This can be written down from symmetry 

and the coefficients and its amplitude fixed by contracting over any two indices). We have also neglected the baryonic pressure 
compared to the radiation pressure. Equations JB3I and ( IB4> form a pair of linear coupled equations for the the perturbations 
in radiation density £q and matter velocity v. Assuming that the rate of variation of of the co-efficients of various terms, due to 
Hubble expansion is small (compared to kc s ), one can use the WKBJ approximation, to derive the dispersion relation for the 
baryon-radiation acoustic oscillations. 

Consider therefore a plane wave solution of the form 



v = V exp(ifc ■ x + Fdr]); 



I exp(ifc • x + / Tdr]) 



(B5) 

Let us also look at longitudinal waves with k parallel to V. Infact, taking the divergence of Eq. ( IB4> one can see that these 
modes are completely decoupled from the rotational modes. To leading order one gets from Eq. JB3I and Eq. (IB4> . a dispersion 
relation which is a cubic equation for V, 
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which can be solved iteratively. Here we have defined b = ( 1 + ?>pb I^Pr) = 1 + R. To the lowest order we get T — ±i(k/ Vhb) . 
So to the zeroth order the dispersion relation is that of a sound (pressure) wave in the baryon-photon fluid, with an effective sound 



speed c s = (1/V36). Consider the effects of terms proportional to L 7 . Since the T 3 term is already multiplied by L 1 we can 



use the lowest order solution to write 



-L 7 r 3 



■L~(— k 2 /36)r. This reduces the cubic equation to the quadratic equation 
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whose solution to first order in L 1 is Eq. ( 16 11 given in the main text. 
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FIG. 2: The figure illustrates the sensitivity of the temperature anisotropy to changes in various parameters, about the standard ACDM model 
of Figure 1. The black (dark) line is always the standard ACDM model. Figure (a): Shows sensitivity to matter/lambda densities for flat 
models. The three plots with decreasing peak amplitudes result from increasing Q, m to 0.5 and 1, keeping total At = 1- Figure (b): Here the 
effect of total curvature on peak location is illustrated. The plot with the smallest I value for the first peak has Q. m — 0.448, Q,a = 0.73, while 
that with the largest is for an open universe, with Q, m ~ 0.27 and Qa = 0. Figure (c): Shows the sensitivity to the baryon density. For the 
three plots with decreasing first peak amplitudes, one has fib = 0.092, 0.046 and 0.023 respectively. Figure (d): The sensitivity to the CDM 
density parameter. This is 0.112, 0.224 and 0.448 for the three plots, where the first peak progressively shifts to lower I values. 
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